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We suggest a long-lived spin polarization structure, a radial spin helix, and study its relaxation 
dynamics. For this purpose, starting with a simple and physically clear consideration of spin trans- 
port, we derive a system of equations for spin polarization density and find its general solution in 
the axially symmetric case. It is demonstrated that the radial spin helix of a certain period relaxes 
slower than homogeneous spin polarization and plain spin helix. Importantly, the spin polarization 
at the center of the radial spin helix stays almost unchanged at short times. At longer times, when 
the initial non-exponential relaxation region ends, the relaxation of the radial spin helix occurs with 
the same time constant as that describing the relaxation of the plain spin helix. 
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I. INTRODUCTION 

At the present time, there is a significant inter- 
est in the field of electron spin relaxation in semicon- 
ductors stimulated by possible future applications of 
spins in electronics and computing^!'!?. In many two- 
dimensional (2D) electron systems the leading mecha- 
nism of spin relaxation is the D'yakonov-Perel' spin re- 
laxation mechanism^;.!. Within this mechanism, electron 
spins feel an effective momentum-dependent magnetic 
field randomized by electron scattering events resulting 
in relaxation of electron spin polarization. A number 
of theoretical and experimental studies on peculiarities 
of D'yakonov-Perel' spin relaxation were reported in the 
last decade^""—. 

It was shown in Ref. 1? that the spin relaxation time for 
2D electrons depends not only on material parameters 
(e.g., strength of spin-orbit interaction, electron mean 
free path, etc.) but also on the initial spin polarization 
profile. In particular, it was demonstrated that a plain 
spin helix in a 2D electron system with Rashba spin- 
orbit interaction has a longer spin relaxation time than a 
homogeneous spin polarization^. A later study ^-''' revealed 
that in a system with both Rashba^^ and Dresselhausi^ 
interactions such an increase in spin relaxation time can 
be even more dramatic. This effect was also observed 
experimentall y ^^1^^ . 

In this paper, we consider spin relaxation of a radial 
spin helix in a 2D electron system with Rashba spin-orbit 
(SO) interaction. This structure is interesting because it 
provides, to the best of our knowledge, the longest spin 
relaxation time of 2D electrons subjected to Rashba SO 
interaction. Such a property is related to spin polariza- 
tion in the vicinity of the special point of radial spin helix 
r — Q. Physically, at short times, when the period of ra- 
dial spin helix is equal to the period of spin precession, 
an electron diffusing from any direction to a point in the 
vicinity of r = has the same direction of spin polariza- 
tion as the initial spin polarization at this point. Thus, 



D'yakonov-Perel' spin relaxation becomes inefficient at 
short times for spin polarization in the vicinity of r = 0. 
In the radial spin helix, the initial distribution of spin 
polarization has a cylindrical (axial) symmetry and is 
given by 



Sr{r,t 



0) = -S'osin(fcr), 
5^(r,i = 0) = 0, 
5^(r,i = 0) = S'ocos(fcr), 



(1) 
(2) 
(3) 



where k is the wave vector, and S'o is the initial amplitude 
of spin polarization. In Fig. [IJa) we show schematically 
spin polarization distribution in the radial direction. The 
overall distribution of z-component of the spin polariza- 
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FIG. 1. (Color online) (a) Schematic of initial spin polariza- 
tion distribution in a radial spin helix, (b) Initial distribution 
of 2-component of spin polarization [Sz (r, Q)/ So) in the radial 
spin helix of finite radius used in our Monte Carlo simulations. 



tion in radial spin helix of finite radius is shown in Fig. 
[IJb). Moreover, in this paper we will often refer to a 
plain spin helix suggested in Ref. M- The initial distri- 
bution of spin polarization components in the plain spin 
helix is 



S^x.t^O) 



-So sin(fca;), 

i}y(X,t = 0) = 0, 

Sz{x, t = 0) = Sq coa{kx). 



(4) 
(5) 
(6) 



This paper is organized as follows. In Sec. [TTl a drift- 
diffusion equation approach is used to study relaxation 
dynamics of the radial spin helix. We obtain general ex- 
pressions for spin polarization as a function of time and 
study its short and long time behavior. Our analytical 
studies are supported by Monte Carlo simulations pre- 
sented in Sec. IIIII Our main results and conclusions 
are summarized in Sec. IIVI Moreover, in several Ap- 
pendices following the main text, we provide additional 
calculation details. Specifically, electron spin rotations 
induced by Rashba SO interaction are considered in the 
Appendix [^ spin drift-diffusion equations are derived 
in the Appendix [B] a general analytical solution of spin 
drift-diffusion equation in the axially symmetric case is 
presented in the Appendix [C] and relaxation of plain spin 
helix is discussed in the Appendix |D] 



II. DRIFT-DIFFUSION DESCRIPTION OF 
RADIAL SPIN HELIX 



about the axis determined by the unit vector n = p x z/p 
(see Appendix [X]) . In this way, the spin-orbit coupling 
constant a enters into equations through the parameter 
1] — 2amh~^ , which gives the spin precession angle per 
unit length. 

Besides this evolution, 2D electrons experience dif- 
ferent bulk scattering events such as, for example, due 
to phonons or impurities. These scatterings random- 
ize the electron trajectories. Correspondingly, the di- 
rection of spin rotation becomes fluctuating what causes 
average spin relaxation (dephasing). This is the famous 
D'yakonov-Perel' spin relaxation mechanism,'^''* The time 
scale of the bulk scattering events can then be character- 
ized by a single rate parameter, the momentum relax- 
ation time T. It is connected to the mean free path by 
(, — VT, where v = p/m is the mean electron velocity. 
To take into account these scatterings we use a model 
of diffusive spin transport, which in the limit of small 
k£ <^1, yields the spin drift-diffusion equations (jB3IIB5[) 
f Appendix [Bl provides derivation details). 

Let us consider dynamics of a radial spin helix relax- 
ation. We assume that such a structure is created at the 
initial moment of time with the spin polarization com- 
ponents given by Eqs. ^^. The exact solution of the 
radial spin drift-diffusion equations (jCHIC2[) with initial 
conditions ([T][3]) and constants 7 and C from Eqs. (|B6p 
can be written as (see Appendix [C] for more details) 



In this section we consider the relaxation dynamics of 
the radial spin helix analytically. The initial spin polar- 
ization in the radial spin helix is of cylindrical symmetry 
and described by Eqs. ^H^. Intuitively, a special point 
in the radial spin helix is r = 0, because the electrons 
motion through this point along straight trajectories in 
all directions should not lead to spin relaxation at short 
times for a specific value of the wave vector k. Corre- 
spondingly, the spin lifetime of electrons located in a re- 
gion within r = should be longer than the spin lifetime 
of homogeneous spin polarization and of plain spin helix. 
This effect is in the focus of our investigation. 

Let us consider a two-dimensional electrons confined 
in a quantum well or heterostructure with Rashba-type 
spin-orbit interaction^^. The standard Hamiltonian with 
the Rashba term is given by 



p2 



iJ = -— + a (5- X p) • z, 
Zm 



(7) 



where p = {px,Py) is the 2D electron momentum oper- 
ator, m is the effective electron's mass, a is the Pauli- 
matrix vector, a is the spin-orbit coupling constant and 
z is a unit vector perpendicular to the confinement plane. 
It is not difficult to show that in the case of Hamil- 
tonian ([7]) the quantum mechanical evolution of a spin 
of an electron with a momentum p can be reduced to 
a spin rotation with the angular velocity Q, = 2ap/h 



Srir^t) _ d r'' dsJiisr) 
So dk Jo 



Vfc2 - s2 



k cosh ( \/rf + l6s^ j 



+ [h] + 4s2) 



sinh fy/rj2 + l6s^nDt\ 



2 J 



So 



dssJo{sr) 



-{s^+37f/2)Dt 



, (8) 



+ (4fc - 7?) 



dkJo Vk^ - s2 



sinh Urj^ + 16s^nDt\ 



cosh I \/rf + l6s^ j 



\/i]'^ + I6s2 



,-{s-+3r,-/2)Dt 



, (9) 



where Ji (r) and Jq (r) are the Bessel functions of the first 
and zeroth order correspondingly and D is the diffusion 
constant. 

Eqs. ([5][S]) define completely the radial spin helix at 
any point r and at any moment of time t. The time 
dependence of spin polarization at the center of helix is 
of particular interest because the spin relaxation at this 
point is the slowest. The radial component of spin po- 
larization Sr in the vicinity of r = is close to zero (it 
follows from symmetry considerations or directly from 
Eq. (HI))- Therefore, below, we derive asymptotic ex- 
pressions at short and long times for Sz only. At short 
times Djft <C 1, an expansion of the RHS of Eq. © in 



t and its integration over s at r = results in 

M^ = l-2ik- TjfDt + -{k - 77)(2fc3 - 6fc2^ 
Do 3 

+6^772 _ 3r]^)D'^f - i—k^ - —k^ 
V 15 5 



^V-^fcV 



+ §k'v'-Y''v' + lv')Dh' + 0{t% (10) 

It follows from Eq. p^ that when the radial spin helix 
period is equal to spin precession length (this happens 
when k = rj) the decay of Sz at the center of helix starts 
with a cubic term in t 



SO 

Co' 



Sz{0,t) 
So 



1 



10, 
"9" 



(Drfty + 0(r), for fc = 77, Dr]H < 1. 

(11) 
This means that the spin relaxation at the center of the 

radial spin helix at short times is significantly suppressed 
(for this special wave number, k — rj) and characterized 
by a rather long initial interval of non-exponential be- 
haviour. 

The asymptotic behaviour of Sr and Sz at long times, 
Drj^t ^ 1, can be determined by taking into account the 
dominant contribution to the integrals in Eqs. ([8|) and 
([9]) . This contribution comes from the vicinity of a point 
s G [0,k] corresponding to the maximum of the — A^(s) in 
the integration interval (see Eq. (|C10p and Fig. [7]). We 
should consider three cases. In the first case, when < 
k < fc„j — \/l5ri/4, the main contribution to the integrals 
in Eqs. ^ and Q comes from the vicinity of point s = k 
at the right end of the integration interval. In the second 
case, when k = fc^, we should keep in mind that k = 
km is a stationary point of A_(s) entering the exponent. 
Therefore, in this case, the asymptotic behaviour differs 
from the case k < km by a pre-exponential factor. In the 
third case, when k > km, the main contribution to the 
integrals in Eqs. ([8]) and ^ arises from the vicinity of 
the inner stationary point s = km of A_(s). 

The asymptotic behaviour of Laplace-type integrals is 
obtained using standard for this purpose technics. From 
Eq. © we get 



Szir,t) 
So 



1 



xe-^-C^'*, for 0< A: 



Ak - T] 
v/772 -f 16A;2 

15 



< 



^ ,, , , Trfc dX_(k) 

d\-(k) , , 

^, -fc— -^t»l,12 
ak 



Sz{r,t) ^15(^15 + 3) /3\ vn2,Ni 

^r ^ 32V2 ^ \V '"^''-'^^^^ *)^ 

xg-A_(fe„)t^ for k = km = ^V, Dift > 1, (13) 



Sz{r,t) 3 (4fc + 5ry)r/ 



So 



32 (fc2 - A:2^)3/2 




Jo{kmr)\I^T-e 



-\-(k^)t 



15 



for k > ^T], {k - kmfDt > 1, Dyft > 1. (14) 
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FIG. 2. (Color online) Time- dependence of Sz(0,t) in the 
radial spin helix at several values of k. 



From Eq. (|T2|) we see that if the wave vector k satis- 
fies the inequality < fc < km, then, at long times, the 
spin polarization decay is mainly exponential and char- 
acterized by a relaxation time T{k) — (A_(A;))^^. This 
relaxation time increases with fc. At the same time (ac- 
cordingly to Eq. (fTU)) ). the spin polarization decay at 
short times also decreases with fc. Therefore, the spin 
life time of the radial spin helix increases with increase 
of fce (0,fc™). 

At fc = km, as it follows from Eqs. (|13m4p . 
the relaxation time reaches its maximum value Tm = 
(A_(fc„))"i = {IDTf/l&y^. Moreover, since fc,„ « 
0.9777 is very close to 77, the conditions for short time 
suppression of spin relaxation are almost optimal at this 
value of fc. Therefore, when fc = km, the spin relaxation 
is significantly suppressed at both short and long times. 

When fc > km, the asymptotic relaxation time is the 
same as when k = km (see Eqs. (|13I14[) '). However, even 
when fc is close to km, the ratio of the absolute value of 
RHS of Eq. HH) to those of Eq. ^ is small (of the order 
of [(fc — km)'^Dt]^^^^). In addition, in the asymptotic 
formula p4|) . the prc-exponcntial factor changes its sign. 
Therefore, the spin polarization Sz must turn to zero at 
some moment of time, before it reaches the asymptotic 
behaviour given by Eq. (|14p . Eq. (fTU]) at fc > 77 w fc^ 
also predicts a relaxation increase with fc when fc > km- 

Thus we conclude that the longest spin relaxation time 
for the spin polarization at the center of the radial spin 
helix occurrs at the wave vector 



1577/4 



and is given by 



{7Dri^/16y 



(15) 



(16) 



In addition, at this value of fc, the dynamics of spin po- 
larization in the vicinity of r = is non-exponential at 
short times, when the spin polarization remains almost 
constant. These are the main results of our calculations. 
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FIG. 3. Spatial dependence of Sz (a) and Sr (b) in the radial spin helix at the indicated moment of time. 



The relaxation of initially homogeneous spin polariza- 
tion (when fc = 0) can be obtained from Eq. ^ in the 
limit fc — > 0. In this limiting case the factor before the 
exponential function g-^-i^)* in Eq. (j9|) turns to zero in 
fc = limit. Therefore, at fc = 0, the spin relaxation is 



determined by the exponential function e 



^A+(0)t 



where, 



accordingly to Eq. (ICTO)) . A+(0) = 2Dtj^. 

In fact, the exact time dependence of 5^ at fc = 
coincides with its asymptotic behaviour. It can be clearly 
seen from both Eqs. © and (jBSP that 



S^it) = 5oe-^+(°)* = 506-2^"'*, for fc = 0. 



(17) 



Accordingly to Eq. (|17l) . the initially homogeneous spin 
polarization, directed along z-axis, decays exponentially 
with a time constant t'''' = {2Dr]^)~^ . We also note that 
the applicability limits of the long times asymptotic ex- 
pressions listed in Eqs. (|12m4l) do not allow calculations 
of 5*2 at fc = or fc = fcm as a limiting case of Eq. (|12|) . 

In Fig. [2] we show the time dependence of the spin 
polarization component Szir = 0,i) calculated at several 
wave vectors (fc = 0.5ry (solid line 1), fc = fc™ = y/T5ri/4: 
(solid line 2) and fc = 1.6ry (solid line 3)) using Eq. ([9|). 
The dashed lines 1, 2 and 3 represent asymptotic behav- 
ior at short times (given by Eq. ([TU])), and dotted lines 
1, 2, 3 show asymptotic behavior at long times (given by 
Eqs. (|12m4p ') for the same values of the wave vectors fc. 
These curves reveal main features discussed above. 

The spatial dependences of z- and r-components of 
spin polarization, calculated from Eqs. ^ and ((H]), are 
presented in Fig. [3] for a particular moment of time 
Drj^t = 10 and wave vectors fc = O.677 (dashed lines), 
k = km (solid lines), and fc — lAij (dotted lines). These 
plots clearly demonstrate that the maxima of Sz and Sr 
are reached at the wave vector fc = km- 

Fig. m depicts time dependence of Sz{r = 0,t) for ho- 
mogeneous spin polarization (calculated from Eq. (jlTI) ). 
plain spin helix (calculated from Eq. (ID4p ) and radial 
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FIG. 4. (Color online) Time dependence of Sz(r = 0, i) for 
three different initial spin configurations: homogeneous po- 
larization, plain spin helix and radial spin helix. The radial 
spin helix is characterized by the largest magnitude of spin 
polarization at any moment of time. 



spin helix. This plot demonstrates that the spin polar- 
ization at r = in the radial spin helix lives longer that 
those in the case of homogeneous spin polarization and 
plain spin helix. It is interesting and important that at 
short times this curve stays almost flat as expected. At 
longer times, when the initial non-exponential relaxation 
region ends, the relaxation of spin polarization in the 
radial spin helix occurs with the same time constant as 
those of the plain spin helix (it can be shown analytically 
from Eqs. (|D3IID4[) that the longest spin relaxation time 
of the plain spin helix is given by Eq. (|16|) and occurs 
at fc = km given by Eq. P^ ). Therefore, in both cases, 
the increase of the exponential relaxation time relative 
to the homogeneous spin polarization is equal to 32/7. 
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FIG. 5. (Color online) Distributions of z component (a) 
and all components (b) of spin polarization [S^p = 0, S = 

{Sl + Siy^) at f = lOOr in the radial spin helix. The red 
(dark) area in the center of radial spin helix in (a) and the 
maximum of S and Sz in the vicinity of r = in (b) demon- 
strate a longer spin lifetime of electrons located in this region. 
This plot was obtained at ?;£ = 0.1 and the radial spin helix 
period a — 64.77£ (this value of a corresponds to fc = km)- 
M.f.p. (mean free path) stands for £. 



III. MONTE CARLO SIMULATIONS 

In order to obtain an additional insight on spin relax- 
ation of the radial spin helix, we perform Monte Carlo 
simulations employing an approach described in Refs. 5 
and [20. This Monte Carlo simulation method uses a 
semiclassical description of electron space motion and 
quantum-mechanical description of spin dynamics (the 
later is based on the Hamiltonian ([7])). All specific details 
of the Monte Carlo simulations program can be found in 
the references cited above and vifill not be repeated here. 
To some extent, Monte Carlo simulations program nu- 
merically solves Eqs. (|B3IIB5I) taking into account rela- 
tions dSl) and (|M|) . 

All numerical results related to the radial spin helix 
were obtained using an ensemble of 10® electrons initially 
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FIG. 6. (Color online) Spin polarization as a function of time 
for different initial spin polarization configurations. This plot 
was obtained using the parameters values rj£ = 0.1 and a = 
64.77^ (the period of helices). The straight lines are fitting 
curves selected as exp( — (f — to)/ti). The parameters of the 
fitting curves are to = 0, ti = 103r for homogeneous spin 
polarization, to = Or, ii = 435r for plain spin helix, and 
to — 165r, fi — 430r for radial spin helix. 



homogeneously distributed within a circle of a sufficiently 
large radius R — 200£ to insure that the influence of 
boundary effects on spin polarization in a region in the 
vicinity of r = is negligible. The initial configuration 
of z-component of spin polarization of these electrons is 
presented in Fig. [Hb). Our Monte Carlo simulations 
results are in an excellent agreement with the theory of 
radial spin helix relaxation presented above. 

Fig. [5fa) demonstrates the spatial dependence of Sz in 
the radial spin helix at t = lOOr and Fig. [5jb) depicts the 
radial dependence of the spin polarization components at 
the same moment of time. In particular, it can be clearly 
seen that the decay of Sz in the vicinity of r = is 
slower than in other regions. The total spin polarization 
S shows oscillations that are related to a well-known fea- 
ture of D'yakonov-Perel' relaxation: the spin relaxation 
time of the perpendicular to plane spin polarization is 
shorter than the spin relaxation time of the in-plane spin 
polarization. Similar oscillation of spin polarization am- 
plitude were previously found in the relaxation dynamics 
of the plain spin helix. ^ 

In addition to the radial spin helix relaxation, we sim- 
ulated the relaxation of homogeneous spin polarization 
and relaxation of plain spin helix. Monte Carlo simu- 
lations reveal that at long times the time dependence of 
spin polarization in all spin configurations (homogeneous 
spin polarization, plain spin helix and radial spin helix) 
exhibits an exponential decay. 

In Fig. [HI we compare the time dependencies of spin 
polarization {S{0,t)/ Sq) for three different initial polar- 
ization configurations: the homogeneous spin polariza- 
tion (initial spin polarization is selected as Sx = Sy = 0, 
Sz = Sq), plain spin helix and radial spin helix. The 



selected period of radial and plane spin helices (fc = km) 
corresponds to the longest spin lifetime of these struc- 
tures at a fixed rjl — 0.1. It follows from Fig. |6] that 
the spin polarization in the radial spin helix is the most 
robust against the relaxation. We emphasize that the 
long-time behavior of all curves can be perfectly fitted 
by an exponential law given in the caption of Fig. |6l 
The numerically obtained increase in the spin lifetime of 
the radial spin helix 435/103 ^ 4.2 is very close to the 
theoretically predicted value 32/7 ~ 4.6. We also note 
that a slightly longer spin relaxation time increase (~ 6) 
reported in Ref. l9| for the plain spin helix can be related 
to a large value of r]£ = 0.3 that is beyond the linear spin 
drift-diffusion theory. 



IV. CONCLUSION 



The dynamics of spin relaxation of radial spin helix 
was investigated using spin drift-diffusion equations and 
Monte Carlo simulations. Starting with a clear model 
of diffusive spin transport, we derived spin drift-diffusion 
equations for electron spin polarization in 2D semicon- 
ductor structures with Rashba spin-orbit coupling. The 
general solution of these equations for the axial symmet- 
ric case was found. Based on this solution, we studied 
the evolution of the suggested long-lived spin structure - 
the radial spin helix. It was shown that the relaxation of 
spin polarization in the vicinity of r = in this structure 
demonstrates an unusual long initial non-exponential re- 
laxation behavior followed by an exponential decay. The 
optimal value of the radial spin helix wave vector was 
found and corresponding exponential relaxation time was 
calculated. Qualitatively, the initial non-exponential de- 
cay feature can be explained by existence of an infinite set 
of dephasing-free trajectories propagating through the 
point r — 0. 

In order to additionally check our analytical results, we 
also performed Monte Carlo simulations of the dynamics 
of radial spin helix relaxation using the same Monte Carlo 
simulation technique as those described in Refs. 5 and 
[20 . Fig. [S]shows a representative result of our simulation 
in which it is clearly demonstrated that the polarization 
decay at r = is slowest. Our analytical and Monte 
Carlo simulations results are in a perfect agreement. 

To conclude, the radial spin helix is a new structure 
exhibiting an unusual spin relaxation dynamics and rel- 
atively long lifetime. Its property of slow relaxation dy- 
namics at short times is very interesting. Experimentally, 
the radial spin helix can be created by spin injection from 
a point electrode located at the center of a second ring- 
shape electrode or possibly by a modified spin gratings 
technique^. 



Appendix A: Rotation of a single spin subjected to 
Rashba SO interaction 



Let us consider an electron with a momentum p 
(px,Py). The Hamiltonian ([7]) can be rewritten as 



H 



P 

2™ 



apn-<T, 



(Al) 



where n = p x z/p is the unit vector. 

Since eigenvalues spin projection operator 
s • n = CT • n/2 on any direction are equal to ±1/2, 
the energy levels of (jAip are given by 



E± 



2m 



±ap, 



(A2) 



and it is clear that corresponding spin eigenfunctions are 
the states with spin directed along and opposite to n. 

The evolution operator acting only on spin variables of 
the electron with momentum p is equal to 



f/(t) 



-iHt/h 



exp 



-ip^t 
2mh 



exp 



-iQtn ■ 



(A3) 
where 51 = 2ap/h. 

We note that the evolution operator (|A3p coincides 
with the operator of finite rotation by an angle fit — 
2apt/h about n-axis up to the phase factor e~'^ t/(2mh)^ 
Since electron spin s transforms under rotations as the 
ordinary vector, we obtain 

s(t) ^ cos(m)s(0)+sin(m)nxs(0)+2sin^(f7i/2)n-s(0)n, 

(A4) 
where 



fi — 2ap/h, n = p X z/p. 



(A5) 



Appendix B: Derivation of spin drift-diffusion 
equations 

Let us consider a two-dimensional non-degenerate elec- 
tron gas and use a semiclassical approach to model 
the electron space motion and quantum-mechanical ap- 
proach based on the Hamiltonian ([7]) to describe the elec- 
tron spin dynamics. Moreover, we assume the electri- 
cal neutrality and absence of an external electromagnetic 
field. Within our approach, 2D electrons are character- 
ized by the momentum relaxation time r and the mean 
free path i, so that the average velocity of electrons is 
V = IJT. From elementary gas-kinetic considerations— 
we can write an equation for the change of electron spin 
polarization AS(x, y, i) in a region of dimensions 21 x 21 
with the center at (x, y) during the time interval r: 

(2£)2AS(a:, y, t) = \vt{2I>) {S'(x - 2£, y, t) 

+S'{x + 21, y, t) + S'(a;, y - 2£, t) + S'{x, y + 2£, t)(Bl) 

-AS{x,y,t)}. 



In the right hand side of Eq. (jBip , the first four terms 
are the spin polarization fluxes into the region from four 
sides with length 2£, and the last term is the flux out 
of this region. The prime symbols in Eq. (jBll) denote 
a change of spin polarization because of SO interaction- 
induced spin precession by the angle 2J7t — Aarai/h — 
2'q(. accordingly to Eqs. (jA4|A5p . For example, 

S'(a; - 2f , y, t) = cos(2??£)S(a; - 2(i, y, t) - sm{2r]£)y 
xS{x-2e,y,t) + 2sin\r^e)y ■ S{x ~ 2i,y,t)y,{B2) 

where y is the unit vector along y— axis. 

In order to obtain drift-diffusion equations for spin 
polarization, we substitute expressions for S' into Eq. 
(|B1I) , and expand trigonometrical functions up to quadric 
terms with respect to small 2r]i and S' terms up to 
quadric terms with respect to 2£. The resulting system 
of drift-diffusion equations for spin polarization have a 
form 



dS, 



dt 



DAS, ~ C 



^ DAS:, 

" = DASy 

dS, 



^^-27^., (B3) 
^^-27^.> (B4) 



dx 



where 



and 



C = 2r]D, 7 



dSy 

dy 



:V^D, 



47^., (B5) 



(B6) 



D ^ 



2t' 



(B7) 



Here D is the coefficient of diffusion, C describes spin 
rotations, and 7 is the coefficient describing spin relax- 
ation. 

It is interesting to note that the same drift-diffusion 
equations (jB3IIB6|) can be obtained for the model of 2D 
localized electrons on a lattice^^ in the hopping regime. 
However, in this case, the diffusion coefficient is equal to 
D = £^/(4t), where r is the characteristic hopping time 
and £ is the distance between lattice sites. 



Appendix C: Analytical solution of drift-diffusion 
equations in the axially symmetric case 

In the axially symmetric case (assuming that Sr — 
Sr{r,t), S^ = S,{r,t) and S^ = 0) Eqs. (IB3IIB5|1 can be 
written as 



dSr 
~dt 



= D 

dS, 
dt 



d 
rdr 



dSr 
dr 



Of 



rdr V dr 



fC^-275.,(CI) 

C^^-475,.(C2) 
ror 



Let us find a general solution of Eqs. (|CmC2p for the 
case of an infinite plane. We search a specific solution of 
the above Eqs. in the form 



Sr{r,t) = A{s,t)Ji{sr), 
SAr,t) ^ B{s,t)Jo{sr), 



(C3) 

(C4) 



where Ji{r) and Jo{r) are the Bcssel functions of the 
first and zcroth order correspondingly. Substituting ex- 
pressions (IC3IIC4I) into Eqs. (ICHIC2I) we obtain a system 
of ordinary differential equations for unknown functions 
A{s,t) and B{s,t) of positive parameter s and time t 



dA{s,t) 

di 
dB[s,t) 

Jt 



= -{Ds^ ^2-i)A{s,t)~CsB{s,t), (C5) 
= ~CsA{s, t) - [Ds^ + A-i)B{s, t). (C6) 



The general solution of this system can be presented as 

A{s,t) =C+(s)sCe-^+('*)*+C_(s)sCe-^- ('*)*, (C7) 
C+(s)(7+^ 



2^g-A+(s)t((^g) 






where we denote 



A± (s) -Z3S2 + 37 ±/ 



r 



(C9) 



and C±{s) are arbitrary functions of positive parameter 
s. Using Eqs. (jB6|) . Eq. (jC9|) can be rewritten in a more 
simple form 

A±(s) = -D{2s^ + irf ± ri^/if + I6s2). (CIO) 

The special solutions (jC3IIC4[) (with A(s, t) and B{s, t) 
given by Eqs. (IC7IIC9I) ') of the radial drift-diffusion equa- 
tions (|CmC2p have the following simple meaning. Ac- 
cordingly to Eqs. (|C3IIC4[) . the spatial dependencies of 
the radial and z-components of spin polarization are pro- 
portional to the first and zeroth order Bessel function 
at any moment of time. The parameter s is similar 
to the wave vector k for the plane case. The ampli- 
tudes A{s,t) and B{s,t) determine the time dependence 
of the radial and ^-components of spin polarization. If 
C+{s) ^ 0,C_(s) = ( C+(s) = 0, C-{s) ^ 0), then 
these amplitudes are exponential functions of time with 
the inverse relaxation time A+(s) (A_(s)) as it can be 
seen from Eqs. (|C7IIC8D . 

Whereas the inverse relaxation time A+(s) takes its 
minimum value at s = and monotonically increases 
with parameter s, the inverse relaxation time A_(s) has 
a minimum at 



(Cll) 





^rn 


= v 


/C2 

4L>2 


1' 


this value of 


s. A- 


.(s)l 


s equal 


to 




Am - 


= 37 


C2 

AD 


7^2^ 
C2 



(C12) 



Using relations (jB6l) . we find that the minimum value 
of A_(s) is equal to 



7 



and it occurs at 



-rj. 



(C13) 



(CM) 



Fig. [7] shows A±(s) given by Eq. (ICIOI) as a function of 
s. 

In order to obtain the general solution of Eqs. (jClt ^ 
IC2I) . we should integrate the special solutions (jC3IIC4[) 
over positive parameter s taking into account relations 
(IC7IIC9p . Two arbitrary functions C±{s) entering Eqs. 
(jC7IIC9[) can be found from specified initial conditions 
for the radial and z-coniponents of polarization in the 
form of its Fourier-Bessel transforms 



Sris) 



S.{s) 



+ 00 



drrSrir,0)Ji{sr), (C15) 



+ 00 



drrS,{r,0)Jo{sr). (C16) 



As a result of algebraic transformations, we obtain the 
solution of the initial value problem (jCmC2p in the case 
of the infinite plane 



Sr{r,t) 



+00 



dssGrr{r,t, s)Sris) 





+00 



dssGrzifji, s)Sz{s), 







+00 



Sz{r^t) = I dssGzrir,t,s)Sr{s) 
Jo 

(■+00 

+ / dssGzz{r,t,s)Sz{s), 
Jo 



(C17) 



(C18) 
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FIG. 7. (Color online) s-dependence of inverse relaxation 
times A±(s). The minimum of A_(s) corresponds to the wave 
vector giving the longest relaxation time. 



where the Green functions Grr, Gzr, Grz, Gzz are defined 
as follows 

Grr{r, t, s) = Mrs) cosh (tVT^TC^) 
sinh (t^/j^ + C^sA ' 



Vt^Tc^I 



-{Ds^+3'r)t 



(C19) 



Grzir,t,s) = -CJiirs) 



Gzr{r,t,s) = -CJoirs) 



sinh (t^/^fTCPl 

X5e-(c-^+37)*^ (C20) 
sinh (t^j^ + C^sA 

x.se-(^'^'+3^)*, (C21) 
Gzz{r,t,s) = Mrs) [cosh (tyf^^+C^) 



'1- 



sinh tvr + c^ 



y/jMC^ 



-{Ds^+3'i)t 



(C22) 



Substituting the initial conditions for the radial spin helix 
(PIS]) into Eqs. (|C15IC16I) and performing integration 
we find that the generalized functions Sr{s) and Sz{s), 
which correspond to the initial conditions ([T]) and ([3]), 
act as follows 



+00 ^ d f^ 

dssF{s)Sr{s) = —So-rr I ds 



d r . kF{s) 



dk Jo Vk'^ - s^ 



dssF{s)Sz{s) — So-jT ds 



sF{s) 
dkJo "" V/c2 1 s2 



,(C23) 



,(C24) 



where F{s) is a smooth enough function. 

Substituting Eqs. (|C23IC24I) into the general solution 
(|C17|C18p and taking into account the expressions for 
the Green functions (jC19IIC22|) . we obtain the explicit 
formulae for the solution of the drift-diffusion equations 
(|CmC2p for the initial conditions (IT]|3t : 



Sr{r,t) = -So- 



ds Ji{sr) 
dk Ja Vk^ - s2 



fccosh 



+ (kj + Cs^ 



sinh (i ^72 + C2s2 



Vr + c^ 

d /■'' dssMsr) 
Sz{r,t) ^So — 



dkj„ V/c2 _ s2 



(fcC-7)' 



sinh (t ^72 + C2f 



e~(D.'+^y)\ (C25) 
cosh (t^/fTc^) 

-{Ds-+3^)t^ (C26) 



V7' + C^i 



Appendix D: Relaxation of plain spin helix 

In the case of the plane spin helix defined by the ini- 
tial conditions (j^HS]), the solution of the drift-diffusion 



equations (|B3IIB5I) has the same harmonical spatial de- 
pendence at any moment of time. Therefore, we can seek 
the solution of this system of equations in the form 

S:,{x,t) = A{k,t)s[n{kx),Sy{x,t) =0, (Dl) 
S:,{x, t) = B(fc, t) cos(fcx). (D2) 

Substituting expressions (|DmD2p into the system (IB31 
IB5|) we obtain the same system of ordinary differential 
equations (|C5IIC6I) for the functions A{k,t) and B{k,t) 
as in the axially symmetric case. Eqs. (jDmD2l) to- 
gether with Eqs. (|C7IIC9p determine solutions of the 
drift-diffusion equations (JB3IIB5I) with such a specific har- 
monic spatial dependence. The arbitrary amplitudes C+ 



and C- are calculated from the initial conditions (Eqs. 
(|3][B])). Finally, the solution of drift-diffusion equations 
(jBSIIBSP for the plane spin helix is given by 



S^{x, t) = -5o sin(A:a:) cosh (^^7^ + C^kA 

.-{Dk-+3j)t^ (D3) 



+ {Ck + 7) 



sinh (iVT^ + C^fc^ 



+ {Ck - 7) 



v/7' + C2fc2 

S^{x,t) = Socos{kx) cosh (tV7^ + C^kA 
sinh (iVT^TC^F) 



v/72 + C^k^ 



-{Dk''+3i)t 



(D4) 
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